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Abstract 

Beam-line equipment was upgraded for experiment E08-027 (g2p) in Hall A at Jefferson Lab. Two beam 
position monitors (BPMs) were necessary to measure the beam position and angle at the target. A new 
BPM receiver was designed and built to handle the low beam currents (50-100 nA) used for this experiment. 
Two new super-harps were installed for calibrating the BPMs. In addition to the existing fast raster system, 
a slow raster system was installed. Before and during the experiment, these new devices were tested and 
debugged, and their performance was also evaluated. In order to achieve the required accuracy (1-2 mm 
in position and 1-2 mrad in angle at the target location), the data of the BPMs and harps were carefully 
analyzed, as well as reconstructing the beam position and angle event by event at the target location. The 
calculated beam position will be used in the data analysis to accurately determine the kinematics for each 
event. 
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1. Introduction 

A polarized ammonia {NHs) target was used for the first time in Hall A for the g2p experiment [T]. It 
operated at a low temperature of IK and a strong transverse magnetic field of either 5 T or 2.5 T. A high 
electron beam current would cause significant target polarization drop due to target temperature rising and 
ionization radiation to the target material [2] . To minimize depolarization, the beam current was limited to 
below 100 nA and a raster system was used to spread the beam spot out to a larger area. The transverse 
magnetic field in the target region would cause the beam to be deflected downward when the beam enters 
the target region. To compensate for this, two chicane magnets were placed in front of the target to pre-bend 
the beam upwards. Due to the low beam current and tight space limitations after the chicane magnets, the 
experimental accuracy goals for the position (1-2 mm) and angle (1-2 mrad) at the target were challenging 
to achieve. New beam-line devices and an associated readout electronics system were designed for the g2p 
experiment to accomplish these goals. Design details and the performance of the beam-line devices will 
be described in the following sections along with a discussion of an analysis method determine the beam 
position and direction. 
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Figure 1: Schematic of beamline components for g2p experiment 


2. Beam-line Instrumentation 
2.1. Beam position monitor (BPM) 

The scattering angle of the outgoing lepton in deep inelastic scattering, which is defined with respect 
to the direction of the incident beam, is an important variable for obtaining meaningful physics results. 
Therefore, the position and direction of the beam, after being bent by the chicane magnetic field and spread 
out by the rasters, must be measured precisely. Two BPMs and two harps were installed for relative and 
absolute measurements of beam position and direction near the target, respectively. 

The BPM consists of four open-ended antennas for detecting the beam position; the measurement is 
non-invasive to the beam. The BPM chambers shown in Fig are part of the beam pipe. The four antennas 
are attached to feedthroughs on the interior wall of the pipe at 90° intervals. The BPM chambers are 



(a) BPM design diagram, from JLab instrumentation 
group 



(b) BPM chamber which contains 4 antennas 


Figure 2: BPM chamber 

placed with a 45° rotation (along z) with respect to the global Hall coordinate. The two pairs of antennas 
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are marked as U- and respectively, which are used to determine beam positions in u and v 

directions. When the beam passes through the BPM chamber, each antenna receives an induced signal. 
The BPM front-end receiver collects and sends the signal to the regular Hall A DAQ system and another 
DAQ system designed for parity violation experiments, the HAPPEX system [3]. The new BPM receiver 
was designed by the JLab instrumentation group [4] in order to achieve the required precision at a level of 
0.1 mm with a beam current as low as 50 nA. The regular DAQ system was connected to a 13-bit fastbus 
ADC (Lecroy ADC 1881) with an integration time of 50 ns, which was triggered by a scattered electron 
event. The HAPPEX system was connected to an 18-bit ADC with an integration time of 875 /is, which 
was triggered by a beam helicity signal at 1 kHz. The amplitude. A, recorded in the ADC has the following 
relation with the BPM signal, (j)\ 


( 1 ) 

where g is the gain of the receiver. 

The BPM receiver generates a large time delay for the output signals. The digital filter used in the 
receiver contributes 1/175 s delay time, which was the inverse of the bandwidth setting chosen for the filter. 
There is a ^ 4 /is delay as a result of finite processing times. The BPM cannot provide event by event 
position because of these time delays, due to the 25 kHz fast raster system. 

Because of the space limitation between the second chicane magnet and the target, the two BPMs were 
placed close to each other. One was placed 95.5 cm upstream of the target while the other was placed 69 
cm upstream, making the distance between them only 26.5 cm. The short distance magnified the position 
uncertainty from the BPMs to target. 

2.2. Super harp 

Two super harps were designed and installed in the beam-line, as shown in Eig{^ (label 6a - 1H04 and 
6b - 1H05A), to provide an absolute measurement of the beam position for calibration of the BPMs. The 
new harps were able to work in pulsed beam (1% duty factor) with a current of several gA. A diagram for 
the harp is shown in EigJ^ which consists of three wires with a thickness of 50 /im, a fork and a controller 



chassis. The harp chamber is perpendicular to the beam pipe and connected to the beam pipe as part of the 
vacuum chamber of the beamline. The two harps have different configurations of three wires: vertical(|), 
bank left(\), and bank right(/) for 1H04, and /, |, \ for 1H05. The angle of the / or \ wire is 45° relative to 
the wire dock frame. The wires are arranged in a fork (EigJ^ controlled by a step motor [5] which can be 
moved in and out of the beam-line. The harps must be moved out of the beam-line when production data 
is being taken because they are invasive to the beam. The original position of the wires was surveyed before 
the experiment at a precision level of 0.1 mm. As the motor driver moved the fork through the beam, each 
wire received a signal, which was recorded for further analysis. The signals received from the wire and the 
step-counters from the motor driver were then sent to an amplifier and the DAQ. The amplification and the 
speed of the motor were adjustable for the purpose of optimizing the signals of each scan. Recorded data 
combined with the survey data were used to calculate the absolute beam position. 
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The signal from the | wire {peak\) was used for getting the x position {x^arp) of the beam, and the signals 
from the /, \ wires {peahj and peak\) were used for getting the y position {yharp)- 

Xharp = survey^ - peak\ 

Vharp = ]^[{survey\ - survey j) - {peak\ - peak/)] (2) 


2.3. Raster system 

In order to minimize the depolarization, avoid damage to the target material from radiation, and re¬ 
duce systematic error for the polarization measurement by NMR (The polarization of the target was 

measured by using a NMR coil which was placed inside the target cell [6] , and the non-uniformity of depolar¬ 
ization could reduce the precision of the NMR measurement due to the measurement being an average over 
the target), two raster systems were installed at ^17 m upstream of the target, as shown in Fig{^ (labels 2 
and 3 for fast and slow rasters, respectively). Both the fast and slow rasters consist of two dipole magnets. 
The same triangular waveforms with frequency of 25 kHz were used to drive the magnet coils of the fast 
raster to move the beam in x and y directions, forming a rectangular pattern of about 2 mmx2 mm, as 
shown in Figj^ 



(a) Magnet current waveform for the fast raster in x 
or y channel 



(b) The 2D histogram of magnet 
current signals of fast raster 


Figure 4: Fast raster pattern 


A dual-channel function-generatoij^was used to generate two independent waveforms to drive the magnet 
coils of the slow raster. The waveforms for the x and y directions are: 


X = Axt^!^ sin(ujk)^ 

y = Ay{t + tQY^‘^sin{ujt^(j)), (3) 

where the Ax and Ay are the maximum amplitude, to and (j) are the AM and sin phase difference between 
X and y waveform, respectively. Both of them are sine functions modulated by a function t^/^ in order to 
generate a uniform circular pattern [7], as shown in Figj^ The frequencies of the x and y waveforms kept 
same: uj = 99.412 Hz. In order to cycle the amplitude modulation (AM) function, four piece-wise functions 
are combined together. The first term is t^/^, and the second term is period — and so on for the third 
and fourth terms. The cycled function has the frequency of 30 Hz. 


^agilent 33522A function generator, http://www.home.agilent.com/en/pd-1871286-pn-33522A/function-arbitrary- 
waveform-generator-30-mhz 
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2 cm 




(a) Magnet current waveform for the slow raster 


(b) The 2D histogram of magnet current 
signals of slow raster 


Figure 5: Slow raster pattern 



(a) Raster pattern with to/0 (b) Simulated raster pattern, with to/0 



(c) Manually adjust to to 0 


Figure 6: None-zero to caused slow raster non-uniformity, (a) and (c) are from the data recorded in the 
ADC, (b) is simulated. The color palette shows the uniformity of the raster pattern. 
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The (f) was locked to ^ by the function generator, while the to was manually fixed to 0. Non-zero to could 
cause a non-uniformity pattern, as shown in FigJ^a), which would cause non-uniformity beam distribution. 
A simulation was reproduced the non-uniformity by setting the to to non-zero, as shown in FigJ^b). The 
to was carefully adjusted and minimized before production data taking to avoid the non-uniformity. The 
pattern of the spread beam was relatively uniform after this adjustment during the experiment, as shown 
in Figj^c). 


3. Data analysis 

3.1. Harp scans for measuring absolute beam position 

An example of a harp scan result is shown in Figj^ There are three groups of recorded data for each 




(a) position vs index for harp scan, used for extend- (b) Signal vs position for harp, x axis is position, y 
ing position record axis is the strength of signal, which is the ADC value. 

Figure 7: 1H05A harp scan data 

harp scan, which are “index”, “position”, and “signal”. The index is related to the moving steps of the fork 
during the scan. Each step of the index increases by 0.008-0.07 mm depending on the speed of the motor 
driver [5]. The position is the wire location for each index. The testing results show a good linear relation 
between the position and the index as shown in Figj^a), because the motor speed is uniform. The line is the 
fitted result with pos = a * index + b. According to this linear relation, interpolation or extrapolation can 
be applied when a few data points are missing, in some cases. The strength of signal vs. position is plotted 
in Figj^b). Each peak represents the location when one of the three wires passed through the beam. 

The positions measured by the two harps were used for calibrating the beam positions in the two BPMs. 
When the chicane magnets were on, beam did not pass straight through from the first harp to the second 
harp. BPM calibrations using two harps were only possible when the chicane magnets were off, i.e. in the 
straight-through settings. Since the BPM was calibrated in the local coordinate system, the calibration 
constants were independent from the settings of other instruments. To make sure that the calibration 
constants for the BPMs were still valid during the non-straight-through settings, the settings for the BPM 
receiver were kept the same as in the straight-through settings during production running. 

The scan data from the harps were not reliable when the current of CW beam (100% duty factor) was 
lower than 100 nA due to the low signal-to-noise ratio. The harp scans were taken in pulsed mode at a 
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current of a few /iA, while the BPMs were used for production data taking in CW mode at a beam current 
of 50-100 uA. For a BPM calibration run, a harp scan was done first in pulsed mode, then a DAQ run was 
taken immediately to record the ADC value in CW mode without changing the beam position. The harp 
scan was then taken again in the pulsed mode to double check the beam position. The harp scan data was 
discarded and the scan was taken again if the beam position changed. 

3.2. BPM data analysis and calibration 

The traditional difference-over-sum (A/S) method of calculating the beam position has the non-linearity 
effect at the position far away from the center of the beam pipe [8]. It is necessary to correct the equation 
of A/S since we have a slow raster with a large size of ~2 cm. With the assumption of an infinitely long 
chamber and neglecting the antenna influence on the electric field inside the chamber, the signal from each 
antenna excited by the beam can be calculated via image charge method (Fig||) [g HOI : 

^ _ 2Rpcos{0i - do) ’ 

where (j)i is the signal received in the antenna, and i is and n_, respectively, is a constant 

related to the geometry of the BPM-chamber and the output resistance, I is the beam current, R is the 
radius of the BPM vacuum chamber, p is the radial position of the beam, and Oi — Oq is the angle difference 
between the antenna and the beam in the polar coordinate . 



Figure 8: Signal for each antenna of BPM 


In order to extract the beam position information, and eliminate the dependence on the beam current 
in equation the A/S method is used as follows: 


Du 


4>U+ (l>U- 


( 5 ) 


where U denotes u and v. Substituting equation into equation (j^, it can be rewritten as follows: 

^ _ 4>u+ - 4>u- _ 2 pcos{0 - 6>o) _ 2 U 

^ <Pu+ + (t>u- R 1+^ -Rl+^’ 

where p^ = v?' v‘^. When v?' v‘^ R^, equation ^ is simplified as: 


(6) 


R _ R - (j)u- 
2 ^ 2 (j)u+ + (j^u- 


( 7 ) 


Equation 0 can be used in the simple case when the beam is near the center of the beam pipe. When the 
beam is far from the center, equation 0 is no longer valid. For the g2p experiment, the beam was rastered 
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• • diff/sum 

‘ * corrected 



-15 -10 -5 0 5 10 15 

real x pos(mm) 


(a) Data from experiment (b) Data from bench test 


Figure 9: BPM non-linearity correction, (a) Comparison between the position calculated from A/S equation 
0 (y axis) and the one from correction equation ^ (x axis). Red solid line is a reference line came from 
linear fit of the center points. Data is collected from the experiment, (b) Comparison between the A/S 
equation 0 and the correction equation 0 using the bench test data. The x axis is the real beam position. 
The red triangles are the positions calculated from correction equation Q. The blue circles are the positions 
calculated from A/S equation Q . 


to have a diameter of about 2 cm at the target. From equation © the beam position is calculated as; 


U 


RDu{ 


Dl+Dl + + 


- 1 ). 


(8) 


The correction equation is tested by using the experiment data and the bench test data. FigJ^a) shows 
the comparison between the position calculated from the correction equation ^ and the one from the A/S 
equation 0 . The red solid line is a reference line came from linear fit of the center points. FigJ^b) shows 
the comparison with the real beam position from the bench test data. In this way the method using equation 
^ can correct the non-linearity effect caused by the A/S method. The handling of the BPM information 
which only used for the center beam position (discussed in chapter [3^ also reduced this non-linearity effect. 

The final information recorded in the data-stream was designed to have a linear response with the raw 
signal in the 50-100nA current range. The (pi in equation © can be rewritten as: 


(pi — Cii{Ai — Ai_ped + ^i), (9) 

where Ai and Ai_ped are the recorded ADC value and pedestal value, and and hi are the slope and intercept 
of the relationship between pi and Ai — Ai_ped- Equation Q can be rewritten as: 


{Auj^ — + ^U+) — hu{Au_ — Au__ped + hu-) 

{Au+ — Au^_ped + ^U+) + hu{Au- — Au-_ped + bu-) 


( 10 ) 


where hu = ajj-jau^^ which is related to the ratio of the signals from the and U- antennas and the 
gain settings of the two channels. 
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Combining the equations and the calibration constant hi was obtained by taking the linear fit 
between the ADC values of BPM and the beam current: / oc {Ai — Ai ped + hi). Besides, the linear fit 
used a group of runs which had the same beam position but different beam current. Figure \To\ shows the 



Figure 10: ADC value of BPM raw signal {A — Aped) V.S. beam current 

~ ^i-ped versus the beam current. It shows that the ADC values were linear with beam current in the 
considering current range of 50-100 nA. The intercept from the linear fit of Figj^is the value hi. 

By transporting the position Xharp and yharp ia equation from two harps to the BPM local coordinate 
Uharp and Vharp^ a fit between the BPM data U and the harp data Uharp determined three calibration 
constants cq, ci and C 2 : 


Uharp = t/c = Co + CiU + C 2 U, (11) 

where Uc is the calibrated BPM position. It was converted to Hall coordinate Xc for further transporting 
to the target location. An calibration example is shown in Fig{^ The asterisks and the dots in FigjTTI 



Figure 11: Harp scan data combined with BPM data, the asterisks are the positions from harp, while the 
dots are from BPM. 
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represent Uharp and [/, respectively. 

In order to reduce the noise and improve the resolution during data analysis, a software filter was applied. 
Since the 18 bit ADC was triggered by the helicity signal with a fixed frequency, it could be regarded as 
a sampling ADC. Fig{T^ shows the signal dealt with a 2 Hz low pass filter. Three plots at the bottom of 
Fig(^ (a,b) are the averaged signal used for comparing with the filtered signal. The results show that the 2 
Hz filter and the 0.5 s average are consistent within the required precision. The filter also erases the beam 
displacement caused by the rasters, which is necessary to extract the position of the beam center. 

3.3. Beam position reconstruction at the target 

It is easy to transport the position from the BPMs to the target by using a linear transportation method 
for the straight through setting. For the settings with a transverse magnetic field at the target, the linear 
transportation method cannot be used since the beam is bent near the target. A target magnet field map 
m was generated from the TOSCA model. To test the accuracy of the TOSCA model, the target magnet 
field was measured before the experiment m usi- The position and angle at target were calculated in 
terms of the positions at BPMA and BPMB, together with the magnet field information. Fits were used 
to speed up the calculation. To do this, a full simulation was taken by generating thousands of trajectories 
with different initial positions and angles. The fits were compared with the full simulation and they are 
consistent with negligible difference. FigfTs] shows the trajectories from the simulation. 

The fitted transport functions were only used to transport the beam center position from the two BPMs 
to the target by applying the 2 Hz filter, which filtered out the fast raster and slow raster motion to keep 
only the beam center position. The transported position were expressed as Xcenter- 

3.4‘ Determining the beam position event-by-event 

The readout of the magnet current for the two rasters was connected to a series of ADCs. Two scintillator 
planes in the HRS form a DAQ trigger. This pulse signal triggered the ADC to record the raster magnet 
current for each event. The information from the rasters and the BPMs was combined to provide the beam 
position event-by-event. The position at the target was determined as: 


— ^center T ^ f straster T ^slraster-) 

where Xfstraster and Xsiraster wcrc the position displaced by the fast raster and slow raster, respectively, 
which were converted from the current values of the two raster magnets. The calibration of the conversion 
factors between the magnet current of the rasters and the displaced position will be discussed in the next 
subsection. An example of reconstructed beam position is shown in Figfl^ 

3 . 4 D. Conversion factor for the slow raster 

Two methods were used to calibrate the conversion factor for the slow raster. The first method used the 
calibrated BPM information, i.e., comparing the raster magnet current with the beam shape shown in the 
ADC of the BPMs. Several calibrations were taken during different run periods at a beam current of lOOnA 
using different values of the raster magnet current, as shown in Fig|T^a). 

The range of the beam distribution at the target was calculated from the ranges at the two BPMs without 
applying the filter, using the transport functions fitted previously. The range of the beam distribution at 
the two BPMs and the amplitude of the raster current were calculated from elliptical fits. An example of 
the fit is shown in Figj^b). Figure [^c) shows a linear fit between the raster current and the range of the 
beam distribution at the target. The x axis in Figjl^c) is the magnet current of the raster, and the y axis 
is the range of the beam distribution obtained from the BPMs. 

The second method for calibrating the conversion factor used a target called “carbon hole” as shown 
in Figl^a). Scattered electrons were used as the trigger for recording the raster magnet current. Since 
the density of the target frame was much higher than that of the “hole”, which was submerged in liquid 
helium, the density of events triggered from the target frame was much higher than that of the hole itself. 
Recorded values reveal a hole shape as shown in Figjl^b). The size of the carbon hole was surveyed before 
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(b) Pedestal run without beam 

11 

Figure 12: Software filter applied to BPM signal, (a) is the signal with beam, (b) is the pedestal signal 
without beam. (1,2,3) in (a,b) are the raw signal without applying the filter, (4,5,6) are applied a 2 Hz 
finite-impulse-response filter with 4th order. (7,8,9) are averaged with 0.5 s. (1,4,7) are the 1-D histogram 
of the recorded signal, (2,5,8) are the signal in time domain, (3,6,9) are in frequency domain. Note all of 
the plots in (a) are from a single signal, same as in (b). 

























































































































































































































Figure 13: Transporting beam position from BPM to target with transverse target magnet field. Trajectories 
are from simulation. Blue lines show the z positions of BPM A, BPMB and target, y and z are in global 
Hall coordinate. 
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Figure 14: Reconstructed beam position at the target 
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(a) ADC value of slow raster, with the raster 
current changing 



(b) Elliptical fit for the spread of mag¬ 
net current of slow raster 



(c) Linear fit between the raster current 
and the range of beam distribution 


Figure 15: Converting the raster current to beam position shift 
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(b) The shape of carbon hole in raster ADC, x and y axis are corresponding 
to the currents on x magnet and y magnet of slow raster, respectively. 


Figure 16: Carbon hole method to calibrate raster 
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the experiment, and a fit program was used to extract the radius of the recorded hole shape for that raster 
current. The conversion factor F was then calculated as the ratio of the size of the carbon hole Shoie and 
the radius of the hole shape Rhoie in the ADC: 


F = 


^hole 
2 * Rhoie 


(13) 


3.4-2. Conversion factor for the fast raster 

The conversion for the fast raster was the same as for the slow raster. The low pass filter for the BPM 
was set to a higher value than the frequency of the fast raster to see the beam shape at the BPM formed 
by the fast raster. For a higher frequency filter, a larger beam current was needed to get a clear pattern. 
The beam current chosen for calibrating the fast raster was near 300 nA, which was the safety limit for the 
target. The beam shape formed by the fast raster is shown in FigfT7| 



11 11 11 11 I 11 11 11 11 11 11 11 11 I 11 11 11 11 11 11 11 I 

-2 -1.5 -1 -0.5 0 0.5 1 1.5 2 

X 


Figure 17: Beam shape formed by the fast raster at the BPM A location, the unit is millimeter 


4. Uncertainty 

The uncertainty of the final beam position at the target for each event contains several contributions: 

• The first part comes from the uncertainty of the calibration constant. It includes the BPM resolution 
for the DAQ runs used for the calibration, the uncertainty of the harp data corresponding to each 
calibration, and the survey uncertainties for the BPMs and harps. It contributes about 0.7 mm for 
the uncertainty of the position and 0.7 mrad for the uncertainty of the angle. 

• The uncertainty on the pedestal is the largest uncertainty for the beam position measurement, con¬ 
tributing about 0.7^ 1.5 mm to the uncertainty of the position and 0.7^ 1.5 mrad to the uncertainty 
of the angle. 

• The uncertainties from the BPM survey need to be included, since the production data and the 
calibration data were taken at different beamline settings when the equipment was moved. They 
contribute 0.5 mm to the uncertainty of the position. 

• The uncertainty from the magnetic field map of the target was considered for the settings with the 
target magnet field. 
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• The uncertainties due to the size conversion of the rasters were also included. 


The position uncertainty was magnified by a factor of 5 at the target because of the short distance between 
the two BPMs. For example, in the straight through setting, if the uncertainty at BPM A is 0.2 mm, and 
at BPM B is 0.27 mm, the uncertainty at the target is 1.1 mm for position and 1.3 mrad for angle. The 
uncertainty for the position at the target was around 1^2 mm, while the uncertainty for the angle was 1^2 
mrad. 


5. Summary 

JLab g2p experiment used a transversely polarized NH^ target for the first time in Hall A. It put a limit 
of below 100 nA on the electron beam current and required a slow raster to spread beam to a large area. 
Two chicane magnets were used to compensate the strong transverse magnetic field. Beam-line equipment, 
including the BPMs, harps and associated readout system, were upgraded to allow precision measurements 
of the beam position at low current (50-100 nA). A software filter was used to reduce noise of the BPMs. 
A correction equation was used to compensate the non-linearity caused by the A/S equation. The harp 
data and the linear fit between the bpm signal and the beam current were used to extract the calibration 
constant of the BPM. To account for the strong target magnetic field effect, transport functions were fitted 
to transport the beam position from the BPMs to the target. The beam position in the x-y plane and the 
angle at the target location are extracted event-by-event by combining information from the BPMs and the 
signals from the rasters. The performance of the new devices (BPMs, harps and slow rasters) were presented 
along with an analysis of systematic uncertainties. 
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